Investigation of the Jahn- Teller Transition in TiF 3 using Density Functional Theory. 



m 
O 
O 

(N 

O 
D 

Q 



O 



Vasili Perebeinos and Tom Vogt 
Department of Physics, Brookhaven National Laboratory, Upton, New York 1973-5000 

(Dated: February 2, 2008) 

We use first principles density functional theory to calculate electronic and magnetic properties 
of TiF3 using the full potential linearized augmented plane wave method. The LDA approximation 
predicts a fully saturated ferromagnetic metal and finds degenerate energy minima for high and low 
symmetry structures. The experimentally observed Jahn- Teller phase transition at T C =370K can 
not be driven by the electron-phonon interaction alone, which is usually described accurately by 
LDA. Electron correlations beyond LDA are essential to lift the degeneracy of the singly occupied 
Ti t2 g orbital. Although the on-site Coulomb correlations are important, the direction of the t2g- 
level splitting is determined by the dipole-dipole interactions. The LDA+U functional predicts an 
aniferromagnetic insulator with an orbitally ordered ground state. The input parameters U=8.1 eV 
and J=0.9 eV for the Ti 3d orbital were found by varying the total charge on the TiFj? - ion using 
the molecular NRLMOL code. We estimate the Heisenberg exchange constant for spin-1/2 on a 
cubic lattice to be approximately 24 K. The symmetry lowering energy in LDA+U is about 900 K 
per TiF3 formula unit. 

PACS numbers: 71.15.Mb,71.15.-m,71.30.+h,64.60.-i 
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I. INTRODUCTION 

There is an ongoing interest in phase transitions in 
perovskite based materials. Above T c «370 K the triflu- 
oride TiF3 has the cubic framework perovskite structure 
AMX3, with no A cations present. Each Ti is at the cen- 
ter of a corner sharing fluorine octahedra MXg. At low 
temperatures the structure becomes rhombohedral. This 
symmetry lowering can to a first approximation be char- 
acterized by a titling of the rigid MX 6 octahedra about 
the threefold axis. The nominal valance of the titanium 
ion is 3+, with one 3d electron occupying the triply de- 
generate t2 g orbital in the high temperature phase. One 
may expect a Jahn- Teller instability in TiF3. Indeed in 
the distorted structure titanium has a D 3t 2 local environ- 
ment, in which the ti g orbitals are split into an a\ g and a 
doublet e g orbital. The trigonal distortions have been in- 
tensively discussed in the context of other ti g compounds, 
like FeO and LaTi0 3 0. 

We use Density Functional Theory (DFT) to investi- 
gate the cubic to rhombohedral structural phase tran- 
sition which has been established by X-ray diffraction 
[10. Many trifluoride MF 3 compounds (M=A1, Cr, Fe, 
Ga, In, Ti, V) exhibit this structural phase transition 
0. It is believed that the observed transition 
in TiF 3 is of a ferroelastic nature 0. In most of those 
materials except for TiF 3 and MnF 3 there are no partly 
filled d shells. The driving mechanism for the ferroelas- 
tic transitions can be pictured as the formation of the 
dipole moments on fluorine atoms due to the asymmet- 
ric distribution of the 2p electron density in the distorted 
structure. The long range dipole-dipole interaction fa- 
vors a distorted structure with anti-polar arrangement 
of dipoles . The Local Density Approximation (LDA) 
captures both the long-range and the short interactions 
on the same footing and predicts a distorted structure to 
have a minimum energy in A1F 3 . 



In the present studies we choose the TiF 3 system with 
a partly filed ti g shell. The Jahn- Teller energy lowering 
due to the lifting of the t2 g orbital degeneracy in the dis- 
torted structure should add up to the long range dipole 
formation energy. Unexpectedly, we find the total ener- 
gies of the high and low temperature phases to be identi- 
cal within the errors of calculations. The failure of LDA 
to explain this phase transition is due to the fact that 
transition metal d electrons are not adequately described 
by LDA. Although the dipole-dipole long range interac- 
tion favors the distorted structure, however the electron 
kinetic energy of the t^ g one third filled band is increased 
due the diminishing of the hopping integral (t) between 
the neighboring Ti atoms due to the octahedra tilting. 
The two effects mutually cancel each other resulting in 
degenerate minimums of the potential energy surface. 

It was realized some time ago that LDA tends to un- 
derestimate the Coulomb repulsion of electrons occupy- 
ing different orbitals of the same d-shell 0. In particu- 
lar, this leads to equal occupations of different orbitals 
of the same manifold and prevents the stabilization of 
the orbitally ordered solutions. The LDA+U functional, 
however, generates an orbital dependent potential which 
favors solutions with broken orbital degeneracy. LDA+U 
calculations indicate that the low temperature phase has 
a by about 900K lower energy per TiF 3 than the high 
temperature phase. One of the roles of the Coulomb in- 
teraction is to suppresses the electron kinetic energy of 
the partly filled t2 g bands. The bandwidth narrowing in 
the distorted structure effectively reduces the restoring 
force to rotate the octahedra which facilitate stabiliza- 
tion of the distorted structure. The electronic ground 
state is orbitally ordered in the low temperature phase 
with one electron occupying the a g orbital oriented along 
the rhombohedral direction forming a Heisenberg spin 
1/2 lattice coupled antiferromagnetically to its neighbors 
via a superexchange mechanism. 
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In this work we use the molecular NRLMOL code 01 
to calculate U by varying the occupation number of the 
Ti d orbital of the TiFg~ ion. These calculations are sig- 
nificantly less computationally demanding compared to 
the supercell approaches |9l llUll2i ]and since only the first 
neighboring fluorines mostly contribute to the screening 
a good estimate for parameters U and J can be readily 
achieved in the calculations. 

The Hund interaction and the electron kinetic energy 
are realistically described with LDA. The TiF3 com- 
pound is predicted to be a ferromagnetic metal with a 
fully saturated magnetic moment of 1 \xb per formula 
unit. The exchange Hund energy is large and the Stoner 
criterium is satisfied. The electron-phonon interaction 
lifts the on site degeneracy of the ti g orbitals and for a 
sufficiently strong coupling a gap in the spectrum can 
open. We estimate the electron-phonon coupling (w 2 
meV/degree) to be insignificant to open the gap and to 
drive the phase transition in TiF3. Although the on- 
site Coulomb correlations are important, the direction of 
the t2g-level splitting is determined by the long range 
dipole-dipole interaction Q ■ Many successes and failures 
of LSDA in transition metal oxide compounds have been 
summarized, for example, in the review article [l3j |. 



II. METHOD OF CALCULATION 



The purpose of this paper is to elucidate the electronic 
structure of TiF3 and explain the structural phase tran- 
sition using the DFT method. This is done using the full 
potential linearized augmented plane wave (FP-LAPW) 
method [14| with local orbital extensions [lj| in the 
WIEN2k implementation 0. The LDA Pcrdcw Wang 
[T?! exchange-correlation potential was used. Well con- 
verged basis sets and Brillouin zone sampling were em- 
ployed. The crystal structure was reported by Kennedy 
and Vogt [4|. It is cubic at high temperature with a 
unit cell volume 58.8 (A 3 ). A rhombohedrally distorted 
structure of the space group i?-3c can be characterized 
by three parameters: (1) volume 56.5 (A 3 ), (2) the octa- 
hedra tilt 4> = 13°, and (3) a c/a — 1.0315 ratio, which 
measures the distortion along the rhombohedral direc- 
tion. There are two titanium and 6 fluorine atoms in the 
rhombohedral unit cell. The fluorine ions are sitting in 
the 6e sites ((x, — x + 1/2, 1/4), etc.), and S — x — 0.75 is 
the deviation from cubic positions. The octahedra tilt- 
ing angle is related to S by tan(</> = 2^/25. The calcu- 
lations for both the high and low symmetry phases were 
performed using two formula per unit cell and the same 
k-point mesh, with 292 special k-points in the irreducible 
Brillouin zone. A tetrahedron method was used for inte- 
gration over the Brillouin zone. 
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FIG. 1: Calculated total energy of ferromagnetic high-T phase 
(4> = 0, c/a - 
c/a = 1.0315 



= 1 - dashed line) and low-T phase 
solid line) as a function of volume. 
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III. LDA RESULTS 

We first calculate the cubic phase by fixing <j> = and 
c/a = 1 in the rhombohedral structure of the space group 
i?-3c. LDA finds a ferromagnetic ground state solution 
with a fully saturated magnetic moment of 1 /is per Ti. 
The volume optimization is shown in Fig. ^ The equi- 
librium volume Vq = 58.6 A 3 is in excellent agreement 
with the experiment £|. The bulk modulus of 111 GPa 
is extracted by fitting a Murnaghan equation of state. 
The c/a ratio was fixed to the experimental value 1.0315 
at T=10 K and the two parameters are relaxed to find 
the minimum of the low symmetry structure. The energy 
minimization with respect to the volume for the fixed 
c/a=1.0315 and <f> — 10° is shown on Fig. ^and yields 
V m = 56.2 A 3 . 

Further optimization with respect to the tilt angle for a 
fixed minimum volume V m and the same c/a ratio did not 
significantly alter the ground state energy. The optimal 
tilt angle is <j>o = 10.5°. The Ti-F bond is the most rigid 
bond in the structure, and therefore it is convenient to 
plot the energy versus the Ti-F bondlength. A parabola 
fit to the energy variation shown on Fig[3 is excellent 
and yields a spring constant of the single bond K = 13.9 
eV/A. The Ti-F bond length d in i?3c crystal structure 
is 



d 2 = a£(l/12 + S 2 + (c/a) 2 /24) 



(1) 



where = a\[2 is the hexagonal a lattice constant. 
Knowing the spring constant for the bond stretching and 
equilibrium volume Vq the bulk modules can be estimated 
as B = K/iW 1 ^) re 96 GPa. 

The total energy differences of the high- and low-T 
phases is beyond the accuracy of calculations. Despite 
the excellent agreement with the experimental structural 
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FIG. 2: Calculated total energy of ferromagnetic low-T phase 
as a function of the tilt angle and fixed volume (V m ) and 
c/a=1.0315 ratio. The optimal Ti-F distance is 1.947 A and 
the spring constant is 13.9 eV/A 2 . 



parameters, the LDA total energy analysis does not ex- 
plain the high temperature of the observed structural 
phase transition. 

The densities of states are shown on Fig. |3| The ti g 
band lies at the Fermi level and it is one third filled. The 
Hund energy dominates the kinetic energy such that the 
Stoner criterium is satisfied and the material becomes 
magnetic. The band structure of the Ti ta g electrons 
can easily be understood with a nearest-neighbor two- 
center Slater-Koster [2(j model. The orbitals \a/3 > have 
hopping matrix elements with themselves along the a and 
(3 directions with amplitude t=(dd7r). The tight binding 
fit in the high temperature phase with t = 0.25 eV is 
shown on Fig. 21 

In the low temperature phase the Ti-F-Ti bond angle 
is smaller than 180°, which reduces the hopping and nar- 
rows the band width. The onsite orbital degeneracy is 
lifted by the octahedra tilting. In the D^d local Ti site 
symmetry the ti g manifold splits into: 
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FIG. 3: Density of states of the ferromagnetic solution for 
high-T phase. The fa 9 and e g bands are split by the crystal 
field energy ~ 2 eV. The t^g bandwidth is about 2 eV, which 
corresponds to hopping parameter t = 0.25 eV. The Hund 
coupling J = 0.87 eV can be estimated from the relative po- 
sition of the spin up and spin down peak positions. 
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FIG. 4: Calculated band structure along the high symmetry 
line in the high symmetry structure: e 9 bands open circles, 
t2 g bands closed circles. The solid line is a tight binding fit, 
which uses a single Slater-Koster parameter t = (ddn) — 0.25 
eV. Fermi level is at zero energy. 
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The bandwidth in the tilted structure <fi = <fi$ corre- 
sponds to reduced t=0.225 eV and the electron-phonon 
Jahn- Teller energy gap at T point is 6 meV. From the 
simple tight binding model with these parameters the ki- 
netic energy loss due to the band narrowing effect is 18 
meV and it is not even compensated by the lifting of 
the orbital degeneracy. The total LDA energy includes 



the lattice and electronic degrees of freedom and pre- 
dicts degenerate minimums within the errors of calcula- 
tions (Fig. P). The LDA failure to predict the cubic- to- 
rhombohedral phase transition can be rationalized by the 
lack of electron correlation effects in the LDA method. 



IV. PARAMETER U (J) CALCULATION FOR 
LDA+U METHOD 

In extending the LDA method to account for corre- 
lations resulting from on-site interactions Anisimov, Za- 
anen, and Andersen (AZA) Q chose to refine LDA by 
including an orbital-dependent one-electron potential to 
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account explicitly for the important Coulomb repulsions 
not treated adequately in the LDA approach. This was 
accomplished in accordance with Hartree-Fock theory by 
correcting the mean-field contribution of the d—d on-site 
interaction with an intra-atomic correction. This correc- 
tion has been applied in slightly different ways. We use 
the SIC LDA+U functional [\'J as implemented in the 
WIEN2k package. 

In this work we use a single TiFg - ion to calculate pa- 
rameters U and J. In these calculations the screening 
due to the nearest neighboring fluorines is taken into ac- 
count and there are no other Ti atoms for the d electron 
to hop to. The TiFg - ion forms an octahedra with a Ti-F 
bond of 1.93 A with the Ti atom placed in the center. 
The total energy and d orbital chemical potential of the 
TiF 6 can be modelled by the single site Hubbard model: 

E = ^o-e(«T+ r H) + ^-K+ n l) 2 -^(«T +n l) 
dE 

Mt = q — = — £ + U(nf + ra j_) — Jn-j- (3) 

where n-f and nj are occupation numbers of the triply 
degenerate ti g molecular orbital. We use two sets of 
self-consistent calculations to determine the parameters 
U and J. The quadratic total energy and linear chem- 
ical potential fits to the nonmagnetic calculations give 
Ui — U—J/2 — 7.66 eV. From the fully polarized mag- 
netic calculation the quadratic energy and linear /ij/2 fit 

we find [/ 2 e// = U- J = 7.24 eV, such that U = 8.08 
eV and J — 0.84 eV. The Hund coupling parameter J 
should be compared with the spin up and down energy 
bands splitting in ferromagnetic LDA calculations Fig. 
El which is J = 0.87 eV. 



V. LDA+U RESULTS 

The volume optimization of the high symmetry phase 
gives Vq = 61.8 A 3 using two formula per unit cell with 
R— 3c point group symmetry operations and fixed (f> = 
and c/a = 1. The bulk modulus extracted from the 
Murnaghan fit shown on Figgis 110 GPa. The electron 
correlations on the Ti d orbitals reduces the bandwidth 
(4t/U), such that the kinetic energy variation with re- 
spect to the lattice constant is smaller in the LDA+U 
functional resulting in a larger equilibrium lattice con- 
stant compared to the LDA result. 

The low symmetry optimization was first done with 
respect to volume for a fixed experimental tilt angle 
= 12° and the ratio c/a = 1.0316. The optimal vol- 
ume is 58.1 A 3 . Further optimization with respect to 
the tilt angle for the fixed equilibrium volume gives an 
additional energy gain of 130 K per Ti such that the 
distorted structure with the optimal tilt angle of 14° is 
lower in energy by 675 K per Ti. The parabola fit to 
the Ti-F bond length variation shown on Fig. gives 
spring constant K = 15.95 eV/A 2 , The corresponding 
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FIG. 5: LDA+U calculations of the antiferromagnetic high-T 
phase (<f> = 0°, c/a = 1 - dashed line) and low-T phase (<f> = 
12°, c/a = 1.0315 - solid line) as a function of the volume of 
TiF3. Coulomb and exchange parameters were chosen U=8.1 
eV and J=0.9 eV. The optimal volumes of 61.8 A 3 and 58.1 
A 3 agrees well with experimental results indicated by arrows. 
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FIG. 6: LDA+U calculations for the antiferromagnetic low-T 
phase as a function of the tilt angle and fixed volume V m and 
c/a=1.0315 ratio. The optimal angle is 14.1°, which corre- 
sponds to Ti-F distance 1.996 A. The spring constant is 15.95 
eV/A 2 . 

bulk modulus B = K/(W 1 /' i ) » 108 GPa is an in excel- 
lent agreement with Murnaghan fit Fig. [3] 

The density of states for the distorted structure is 
shown on Fig. [7| The fluorine 2p states are not much af- 
fected by the correlations. Whereas dramatic differences 
are seen for the Ti d states. First, LDA+U predicts a 
gap for both high-T and low-T phases in the ti g band, 
which is split into three distinct narrow peaks originated 
from the a g and two e g orbitals. The imposed rhombohe- 
dral symmetry of the supercell makes the lowest energy 
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a g orbital pointing along the (111) direction to be fully 
occupied with a spin up electron on one atom and a spin 
down on the other atom coupled antiferromagnetically. 
The electronic properties of TiF3 can be described by 
the Hubbard Hamiltonian: 



t-f 2 

+ £ 

i,a(3^a' f3' 



U' - JS a . 



(4) 



The first term is a kinetic energy term described in Sec- 
tion [ffi] To place two electrons with opposite spins costs 
energy U, if they occupy the same orbital, or U', if they 
are on different orbitals. If the spins of two electrons on 
different orbitals are aligned the energy cost is reduced 
by J. In the atomic limit the relation U = U' + 2J 
holds. Then the splitting between the a g and e g orbitals 
for a single site Hubburd model would be U' — J » 5.6 
eV. The atomic limit value of splitting is larger from the 
complete LDA + U solution w 4 eV (see Fig. [7J where 
only a fraction of d electrons inside the muffin-tin sphere 
experiences the LDA+U potential. The second order per- 
turbation theory of the Hamiltonian (Eq. with respect 
to t/U predicts an insulating antiferromagnetic ground 
state solution for parameters U = 8.1 eV, J — 0.9 eV 
and t = 0.22 eV. 



The superexchange energies per Ti of the antiferro- 
magnetic Neel state and ferromagnetic solutions, with a g 
orbital Eq. |21 occupied in both cases, are: 



Eafm — 
E FM = 



4t2 

IT 



1 



8t 2 1 

~T~u 



3 U'-J 



(5) 



where the first term is due to virtual hopping of the lo- 
calized electron on the six neighboring e g orbitals. The 
second term in Eafm is due to virtual hopping of the a g 
orbital to the 6 neighbors with the hopping amplitude 
24/3, and it is missing in the ferromagnetic ground state 
due to the Pauli principle. The LDA + U energy differ- 
ences between the ferromagnetic and antiferromagnetic 
ground states and the fixed optimal distorted geometry 
are shown on Fig. |SJ Both the perturbation theory and 
the LDA + U results predict the antiferromagnetic to 
ferromagnetic transition at U / J ~ 4. For the large U/J 
values LDA + U underestimates the exchange energy by 
a factor of 2 and for the realistic parameter of U = 8.1 
eV it predicts an energy difference of 72 K between the 
two phases. Neglecting the zero-point energy of the spin 
waves the energy difference between the parallel and anti- 
parallel classical spins on a cubic lattice is 6JS(S + 1) 
J23j. For spin S = 1/2 this yields J = 16 K and the cor- 
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FIG. 7: Density of States per two formula units per eV. The 
two electrons in the unit cell fill the a} g spin up orbital on 
atom one and the spin down orbital on the second atom 
peaked at -1.64 eV. The two lowest unoccupied peaks (at 2.01 
eV and 2.52 eV) are due to the e g orbitals of the tig manifold. 
The lower energy orbital is due to the e g orbital of the same 
atom and spin as the occupied a g orbital. The broad bands 2 
eV higher in energy shown by the dashed and dotted lines are 



'} orbitals split by crystal 



due to the e g = {x 2 — y 2 ,3z 2 
field. The empty a g orbital with opposite spin at 4.53 eV 
forms a resonance in the crystal field e g background. 



responding Neel temperature in the mean field approxi- 
mation is Tn = 1.5 J = 24 K. The quantum fluctuations 
reduces the Neel temperature to Tn = 0.946 J = 15 K 
[24|. However, the neutron diffraction experiment [2{| 
did not reveal any long range magnetic order down to 
10 K. Whether the absence of magnetic long range or- 
der in TiF3 is an intrinsic effect due to the coupling of 
the orbital and spin degrees of freedom or extrinsic due 
to the presence of the multidomain structure observed in 
the low symmetry phase by Mogus-Milankovic et al. Q 
requires further theoretical and experimental investiga- 
tions. 



VI. PRESSURE ANOMALY 

X-ray diffraction measurements under pressure by 
Sowa and Ahsbahs [26j allow us to determine the ex- 
perimental bulk modulus of TiF3. Fig. El (a) shows data 
points along with Birch-Murnaghan equation of state fit. 
The fit to the full range of pressures predicts a physically 
unreasonable small bulk modulus B = 7.5 GPa and un- 
physically large coefficient K2 — 142 GPa -1 . The first 
derivative is fixed to K\ = 4 in the Birch-Murnaghan 
equation of state: 



P(V) = 3S/(1 + 2/) 5 / 2 [l + -^[BK 2 + ^ 
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FIG. 8: Energy difference of the total energy of ferromagnetic 
and antiferromagnetic solutions as a function of U/J ratio. 
The solid line is a second order perturbation theory for the 
fixed values J=0.9 eV and t=0.225 eV. The dots are LDA + U 
results with fixed J = 0.9 eV and optimal low temperature 
geometry. 
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FIG. 9: (a) experimental P — V diagram fitted by the Birch- 
Murnaghan equation of state to all data points (solid line) 
and the first six high pressure points (dashed line), (b) bulk 
modulus calculations using the cell parameters given in table 
HI Energies of LDA (LDA+U) calculations are shown by open 
(filled) circles. The Birch-Murnaghan fit to the high pressure 
points is shown by the dashed (solid) line for LDA (LDA+U) 
calculations. The zero energy is chosen at the ambient pres- 
sure value. 



where / = 0.5((Vb/V) 2 / 3 — 1) and Vb is the volume at am- 
bient pressure. The high pressure points (above 4 GPa) 
can be fitted by Eq. JBJ to give B = 51 GPa and a 
normal K 2 = —0.005 GPa" 1 . However, the volume at 
normal pressure has to be fixed to Vb — 51.9 A 3 , which 
is 11% smaller than the observed one. 

In the present work we calculate the bulk modulus of 
the lower symmetry phase using the experimental c/a 
and V/Vo ratios [26(. The tilt angle was chosen to keep 
the Ti-F bondlength Eq. constant for all pressure 
values. We fix d T i-F = 1-947 A and d T i-F = 1-996 
A for LDA and LDA+U calculations respectively. For 
LDA+U calculations Vq was fixed to the experimental 
value 26] V =58.3 A 3 , while for LDA it is reduced to 
match the optimal volume Vb =56.2 A 3 to eliminate the 
strain effects. The results of calculations for the unit cell 
parameters given in tableware shown on Fig. El(b). The 
high pressure points can be fitted by Eq. © to give a 
bulk modulus of B = 42 GPa for LDA and B = 29 GPa 
for LDA+U functionals with corresponding coefficients 
K 2 = -0.011 GPa" 1 and K 2 = -0.008 GPa" 1 and op- 
timal volumes V/Vq = 0.86 and V/Vq — 0.88, which are 
in agreement with the ambient pressure volume deduced 
from the high pressure experimental data. The energies 
in the global minimums in LDA and LDA+U calcula- 
tions are 270 K and 525 K per TiF3 lower than that of 
the experimentally observed structural parameters. A 
shallower minimum in LDA can be explained by the ki- 
netic energy loss due to the octahedra tilting. In LDA+U 
the kinetic energy is suppressed due to correlations and 
does not change much with tilting. 

We speculate that this anomalous behavior of the com- 
pressibility of TiF3 may be due to the presence of domain 



TABLE I: The unit cell parameters for the bulk modulus 
calculations. The normalized volume V/Vo and c/a ratios 
are taken from Sowa and Ahsbahs [2fil |. For LDA+U calcula- 
tions we used experimental volume Vb =58.3 A 3 at ambient 
pressure |26| and we choose a smaller Vb =56.2 A 3 for LDA 
calculations to match the LDA optimal value. For given vol- 
ume and c/a ratios the tilt angle was chosen such as Ti-F 
bondlength Eq. Q is constant and equal to dn-F = 1.947 
A and d T i-F = 1.996 A for LDA and LDA+U respectively. 



p (GPa) 


V/Vo 


c/a 


<f>° LDA 


<j>° LDA+U 


0.0001 


1.000 


1.023 


10.4 


13.9 


0.46 


0.957 


1.056 


14.3 


17.0 


1.05 


0.908 


1.099 


17.9 


20.2 


1.62 


0.885 


1.118 


19.4 


21.5 


2.37 


0.864 


1.138 


20.8 


22.8 


3.18 


0.848 


1.155 


21.7 


23.6 


4.03 


0.830 


1.167 


22.8 


24.6 


4.97 


0.818 


1.181 


23.5 


25.3 


6.05 


0.806 


1.187 


24.2 


25.9 


6.40 


0.801 


1.189 


24.4 


26.2 


7.23 


0.793 


1.193 


24.9 


26.6 


7.67 


0.788 


1.198 


25.2 


26.9 



walls in the low symmetry phase observed in [3j. Under 
applied pressure the crystal becomes a single domain. 
An alternative explanation is the phonon contribution to 
the free energy which in principle depends on pressure. 
Additional theoretical calculations of the phonon spec- 
trum and it's pressure depends (Gruneisen parameter) is 
needed to estimate the phonon contribution to the bulk 
modulus anomaly. 
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VII. CONCLUSION 

Wc used density functional theory to predict the elec- 
tronic and magnetic properties of TiF 3 . LDA predicts 
TiF 3 to be a ferromagnetic metal with a fully saturated 
moment 1 /is per Ti. The energies of the high and 
the low symmetry structures are degenerate. Such that 
the pure electron-phonon and electrostatic (Madelung) 
model does not explain the observed phase transition at 
T c =370 K. The correlations are essential to suppress the 
kinetic energy loss of Ti t2 g electrons to favor the dis- 
tortions. To model electron correlations on Ti d orbitals 
we use the LDA+U approach, which requires the input 
parameters U and J. We determine these parameters by 
calculating electron correlations on the TiF|?~ ion and 
find U = 8.1 eV and J = 0.9 eV. LDA+U predicts TiF 3 
to be an antiferromagnetic insulator with spin 1/2 per 
Ti. We find a long range order of Ti a g orbitals with a 
wavevector (000). The low temperature phase is lower 
in energy by about 900 K per TiF 3 in LDA+U calcu- 
lations, which suggests electron-electron correlations are 
important. 

Using the experimentally determined c/a and V/Vq ra- 



tios we find a global minimum at 14% (LDA) and 12% 
(LDA+U) volumes smaller than the observed ambient 
pressure volume, which are consistent with the ambi- 
ent pressure volume deduced from the experimental high 
pressure data. The experimentally observed larger vol- 
ume at the ambient pressure and the unusual behavior of 
the bulk modulus could be understood once the phonon 
contribution to the free energy is included. The zero 
phonon energy is in principle pressure dependent and 
can alter the position of the energy global minimum and 
the bulk modulus. The Grunciscn parameter calculations 
and experimental studies of the phonon spectrum under 
pressure will help to elucidate light to this problem. 
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